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Abstract 



We study the dynamics of a polymer when it is quenched from a 9 solvent 
into a good or bad solvent by means of a Langevin equation. The variation of 
the radius of gyration is studied as a function of time. For the first stage of 
collapse or swelling, the characteristic time-scale is found to be independent 
of the number of monomers. Other scaling laws are derived for the diffu- 
sion regime at larger times. Although the present model is solved only for 
homopolymers and doesn't include hydrodynamic interactions, these results 
may be a first step towards the understanding of the early stages of protein 
folding. 

Under certain conditions (such as a temperature decrease or exposure to a bad solvent), 
homopolymers undergo a collapse (coil-globule) transition. This has been studied quite 
extensively, and it has regained some interest since it may be the simplest model to describe 
the first stages of protein folding. Few experimental data are available JTJ 0, but a lot of 
theoretical work and computer simulations have been done, leading to a large diversity of 
results. Although the thermodynamic properties of the homopolymer problem have been 
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quite well understood, there is at present no analytical derivation of the kinetics laws for the 
polymer collapse. Since one might expect the first stages of protein folding to be dominated 
by non-specific hydrophobic forces, this phenomenon should be quite similar to the collapse 
of a homopolymer chain. It is thus highly desirable to develop a general analytical technique 
in order to study the kinetics of polymer collapse. 

According to de Gennes' theory ||, the process of collapse of a flexible coil leads to the 
formation of crumples on a minimal scale along the linear chain, which thickens and shortens 
under diffusion of the monomers, then forms new crumples of growing scale, until the final 
state of a compact globule is reached; the total time of collapse is estimated as 

r]a 3 \AT\ 
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where rj is the viscosity of the solvent, 9 is the 9 temperature, a is the monomer size and 
AT is the temperature quench. This time r c has a very strong dependence on molecular 
weight. For the case of proteins (see reference 0), N — 300, which yields a collapse time of 

T c ~ 1 /IS. 

In a revised model, A.Buguin et al. [|J consider a different mechanism, in two steps: first 
a fast formation of "pearls" along the chain, followed by a slower stage of compaction. 

Grosberg et al. [f| have included the role of self entanglement to de Gennes' model, 
and consider a two stage mechanism: first a collapse leading to a crumpled globule with 
characteristic time ~ N 2 , and then a chain knotting driven by reptation with a longer 
characteristic time ~ iV 3 . 

Monte Carlo studies (Ostrowsky&Bar-Yam ||, Milchev&Binder JF|) and Langevin dy- 
namics simulations (Byrne et al. ||) lead to diverging interpretations concerning the mech- 
anisms of collapse. 

In a series of articles, Timoshenko et al. || have developed an alternative theory based on 
a self-consistent method using Langevin equations that can be analyzed numerically; kinetics 
laws for the collapse of a homopolymer are obtained with or without hydrodynamics, at early 
and later stages. 



Experimentation in this field is quite difficult since one has to work in a very dilute regime 
in order to avoid aggregation of chains and observe a collapse. Experiments by Chu et al. |I[ 
on polystyrene in cyclohexane reveal a two-stage kinetics of approximately equal times. The 



folding by ultra-rapid mixing; based on optical techniques, these experiments can monitor 
folding up to the microsecond time-scale. 

In the following, we shall present an analytical method to study the kinetics of a ho- 
mopolymer in a 9 solvent when it is quenched into good solvent conditions (swelling into a 
coil) or into bad solvent conditions (collapse into a globule) . 

We consider a homopolymer chain in a 9 solvent - i.e a Gaussian coil- consisting of N 
monomers, obeying the Langevin dynamics as the chain is quenched into good or bad solvent 
conditions (equations (jl|) and @ ). 

Neglecting hydrodynamical interactions (the Oseen tensor) the equations of motion for 
the system read: 



where N is the total number of monomers, r(s,t) is the position of monomer s in the 
chain, a® is the monomer length and Tq = where D is the diffusion constant of a monomer 
in the solvent and ksT is the temperature. The intra- molecular as well as intermolecular 
interactions of the chain are contained in the potential V(r(s,t)). The thermal noise t](s,t) 
is a Gaussian noise with zero mean and correlation given by: 



The method consists in finding a virtual homopolymer chain which obeys a simpler 
Langevin equation, chosen so that its radius of gyration best approaches the radius of gyra- 
tion of the real chain at each time t. 

The virtual chain, defined by r^(s, t) satisfies the Langevin equation: 



most promising experiments are those by Chan et al. |2j who study sub-millisecond protein 




(1) 



(2) 



< rj(s, t)r](s, t') >= 2D8(t - t') 
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with the same diffusion constant and noise as the original equation, but with a much 
simplified Hamiltonian H v . Indeed this Hamiltonian H v represents a Gaussian chain, but 
with a time dependent Kuhn length a(t). 

Our method is a generalization of Edwards' Uniform Expansion Model Jl0| to dynamics. 
This method consists in calculating the radius of gyration of a polymer in a good solvent 
by using perturbation theory, and adjusting the simplified Hamiltonian so that the first 
order perturbation to the radius of gyration vanishes. If v denotes the excluded volume, the 



method gives the Flory radius |TT| for large N and agrees with the result of the first-order 
perturbation expansion for small v. Note that it would seem natural to use the most general 



quadratic Hamiltonian rather than that of (^) , but this was shown by des Cloizeaux |12j to 
yield the incorrect exponent v — 2/d. 
Let's define 



X(s,t)=r(s,t)-r^(s,t) 
W = H -H v 

Assuming that ([3]) is a good approximation to (|1|) , x(s,£) an d W can be regarded as 



small, and to first order in these quantities, the dynamical equations become: 



dt ~ a 2 (t) ds 2 

dx r d 2 x 



dt a 2 (t) ds' 



+ r 



Q2 r {v) 



v Oq a? (t) ) ds 2 



+ F(rW{s,t)) 



(5) 
(6) 



where F(r(s, t)) 



dv 



dr(s,t) 



is the driving force for the swelling or collapse of the chain. 



In this approximation, the radius of gyration of the chain becomes: 



1 r N 

R g = J o < r 2 (s,t) > ds 



1 rN 

-J Q < ((r^) 2 (s,t) + 2r^(s,t) X (s,t)) > ds 



(7) 
(8) 



The brackets denote the thermal average ( that is an average over the Gaussian noise 
r](s,t)). Our approximation consists in choosing the parameter a(t) in such a way that the 
first order in (JS]) vanishes: 

r N , \ 

\ <r^(s,t) X (s,t)>=0 (9) 
J o 

or in Fourier coordinates: 

Y,<r ( : ) m:(t)>=o (io) 



where the Fourier transform is given by: 



We have used periodic boundary conditions, so that u n = In addition, to get rid of the 
center of mass diffusion, we constrain the center of mass of the system to remain at fixed 
position, r (t) = ^'(t) = Xo(t) = 0. 

Equations @ and (]|) can easily be solved in Fourier space. We assume that at time t = 0, 
the chains are in a 9 solvent, so that the initial condition {r(s, 0)} obeys Gaussian statistics. 
We choose the initial virtual chain to coincide with the real one, so that r^(s,0) = r(s,0) 
for any s. Denoting by TTT the average over the initial conditions, the correlation function 
of r(s, 0) (in Fourier space) is taken as: 



f„(0)= (11) 
U0)r*(0)= ^- 2 5 mn (12) 

Replacing r$(t) and x*(t) by their expression in (|T0l), and taking thermal and initial 
condition averages, we obtain an implicit equation for a(t). 

This equation can be solved analytically in both limits t « tr (short time limit) and 
t » tr (long time limit) where tr = -^r^ is the Rouse time. 
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In the following, we neglect all hydro dynamic interactions and consider only two-body 
and three-body interactions between monomers. 

More precisely, for a chain in a good solvent, we consider excluded volume interactions: 



V(r(s, t)) = V 2 (r(s, *)) = -/ ds ds'6(r(s, t) - r(s', t)) 

2 Jo Jo 

and for a chain in a bad solvent, we take attractive two-body interactions and repulsive 
three-body interactions: 

V(r(s,t)) = -V 2 (r(s,t)) +V 3 (r(s,t)) 



V(r(s,t)) = -~ ds dJ6{r(s,t)-r(sr,t)) 
lo Jo 



V 

2 

w 



rN rN fN 

+- ds ds' ds"5{r{s,t)-r{s , ,t))6{r{s',t)-r{s",t)), 
6 Jo Jo Jo 

where v > and w > 0. 

In good solvent conditions, at short times t << Tr, the radius of gyration increases like 
a power law, which we recast in the form of a stretched exponential: 

R 2 Jt) = Na 2 e^, (13) 



where the characteristic time r c is defined below in (|i9|) . 

For large time t » tr, the radius of gyration relaxes to its Flory value: 

R g (t) ~N^a (l -e"^), (14) 

with a relaxation time t\ given by: 

d+S 

T\ ~ N d + 2 

Note that this relaxation time is much larger than the Rouse time, for dimensions lower 
than 4. For example in d = 3, t\ ~ . 

In a bad solvent, at short times t << tr, the radius of gyration decreases as above: 



R 2 (t)=Na 2 e-V-, (15) 



with a characteristic time r c defined in ( |T9|) and for large times t >> r^, the radius of 
gyration relaxes to that of compact globule according to 

R g (t)~(^N*(l-e^), (16) 

where 

,w. 2 , , 2 

r2 rv, (_)djV 1+ d. 

f 

Note that for dimensions larger than 2, the relaxation time is much shorter than the 
Rouse time. For example in d = 3, T2 ~ Na. 

In the two cases discussed above, the dynamic relaxation exponent z, relating the relax- 
ation time of the system to its radius of gyration: 

r ~ W g (17) 

satisfies the exponent relation: 

2 = 2 + ^ (18) 

This relation was derived by de Gennes fl"3|| . 

The first stage of the swelling or collapse can be described by a stretched exponential 
with characteristic time r c given by: 

ri = - — ° (19) 



where 



3 V2D^vI d N^ 



l rl 1 

g?u / rfw j with \u — u'\ > A 

o Jo flu-w'Kl - \u-u'\)]2 



and A = 1/iV is a short distance cut-off. 

For d < 2, the integral converges for small A, and Id is independent of N. On the other 
hand, for d > 2, the integral is infra-red divergent and thus there is an explicit dependence 
on the cut-off. It is easily seen that this N dependence exactly cancels out the dependence 
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in ( |19D so that the final characteristic time r c is finite (independent of N). In particular, for 
d — 3, we find 



2tt 2 /a^ z / - n2 



Tc ~ 9D V u / V5.22 



The typical order of magnitude of this short time collapse can be calculated for a 
typical protein in water. The diffusion constant of an amino-acid in water is typically 
D ~ 10~ 7 cm 2 /s . Taking a monomer length of ~ 4A, for a chain of 100 aminoacids, 
we find a microscopic characteristic time r c ~ 10 _9 s (in these conditions, the Rouse time 
is tr ~ 4.10 -6 s). Note that in these conditions, the relaxation time T2 is of the order of 
magnitude of the Rouse time. This time is several orders of magnitude lower than other 
estimates in the literature ( see references [f[4j| , [fj] ). The fact that this characteristic time 



is independent of N shows that in the early stage, the collapse is a very local phenomenon, 
where nearby monomers aggregate into small domains. 

Whereas previous mechanisms relied mainly on pure diffusion, we have introduced here 
explicit attractive hydrophobic forces for the collapse of the chain. These forces are deter- 
minant at the early stages of the process, and make it much faster than pure diffusion. 

At larger time scales, we find other scaling laws for the relaxation towards equilibrium 
(which is either the collapsed globule or the swollen coil within the Flory theory). Let us 
mention that we find the same relaxation time for the collapse as in Kuznetsov et al. ||, 
namely ( T2 ~ iVa). This time is several orders of magnitude larger than r c 

Although in this work we did not take into account the role of hydrodynamic interactions, 
we expect their effect to be weak, at least at the early stage of the collapse. However, these 
effects might become important in later stages of folding. 
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